' LyZj ^ WlL vovw 



NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 2751 



r 

h 

Tt 


A SIMPLE APPROXIMATE METHOD 

FOR CALCULATING SPANWISE LIFT DISTRIBUTIONS 

AND AERODYNAMIC INFLUENCE COEFFICIENTS 

AT SUBSONIC SPEEDS 

By Franklin W. Diederich 

Langley Aeronautical Laboratory 
Langley Field, Va. 



Washington 
August 1952 


— 

TECHNICAL ! 'BCARY 
AFL 2811 


TECH LIBRARY KAFB, NM 




NATIONAL ADVISORY COMMITTEE FOR AERONAUT 


TECH LIBRARY KAFB, NM 


illllilil 

DDbSBfl 


TECHNICAL NOTE 2751 


A SIMPLE APPROXIMATE METHOD 
FOR CALCULATING SPANWISE LIFT DISTRIBUTIONS 
AND AERODYNAMIC INFLUENCE COEFFICIENTS 
AT SUBSONIC SPEEDS 
By Franklin W. Diederich 

SUMMARY 


Several approximate methods for calculating lift distributions at 
subsonic speeds are combined and extended to form a simple step-by-step 
procedure for calculating symmetric and antisymmetric lift distributions 
for arbitrary angle-of-attack conditions on swept and unswept wings. 
Methods of estimating the required aerodynamic characteristics are 
included, but any available theoretical or experimental results may be 
used in several steps of the analysis to shorten the work and increase 
the accuracy. The extension of the method to the calculation of aero- 
dynamic influence coefficients and of spanwise moment distributions is 
indicated. 


INTRODUCTION 


An empirical method (ref. l) for calculating spanwise lift distri- 
butions on unswept wings has been used extensively in the past. In 
modified form (ref. 2 ) it has been applied to the calculation of span- 
wise lift distributions on swept wings. More recently, an improved 
method of calculating the lift distributions due to twist has been 
published (ref. 3 ). 

In the present paper a limitation of the method of reference 3 as 
applied to antisymmetric twists is indicated and removed. The methods 
of references 1 to 3 are then combined with each other and with the 
results of another investigation (ref. 4) into a step-by-step procedure 
which starts with the lift-curve slope and the additional lift distri- 
bution and proceeds with the calculation of symmetric basic lift distri- 
butions, the lift distribution due to roll, the rolling -moment coeffi- 
cient due to roll, and, finally, antisymmetric lift distributions. The 
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additional lift- distribution is herein defined as the lift distribution 
for constant angle of attack across the span with a lift coefficient 
equal to 1. The basic lift distribution is defined as the lift distri- 
bution of a twisted wing airzero (total) lift. Means of calculating 
the aerodynamic characteristics required in the method are contained in 
this paper, but if any of these characteristics are known from other 
sources they may be incorporated in the procedure with a resultant 
saving in t im e and improvement of accuracy. 

The method of this paper is derived in appendix A. It is outlined 
and some examples of its application to various plan forms and angle-of^ 
attack conditions are presented and discussed in the body of the paper. 
Formulation of the method in matrix notation for the purpose of obtaining 
aerodynamic influence coefficients suitable for aeroelastic analyses is 
accomplished in appendix B of the present paper in the manner employed 
in reference 5- Inasmuch as the method of the present paper supersedes 
that of reference 2, on which reference 5 is based, appendix B of=the 
present paper supersedes reference 5* 

Means for estimating the spanwise distributions of pitching moments 
(or local centers of pressure) required in an aeroelastic analysis are 
indicated in appendix C. 
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SYMBOLS 

aspect ratio, b^/s 

distance of local aerodynamic center from leading edge, 
fraction of chord 

coefficient in Fourier series for 7 

angle of attack, radians unless specified otherwise 

average angle of attack, radians 

effective angle of attack, radians 

induced angle of attack, radians 

dc 7 /dc 7 

control-effectiveness parameter, •/ 

dS / da 


b 


wing span 
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*n 

C l,2,3 

c l 

Cz 

c *d 

c h 

c 

c 

c l 



C P8 


7 

?b 

?d 

?r 

?1 


coefficient in Fourier series for a sin © 
coefficients for additional lift distribution 

wing lift coefficient, Lift/qS 

wing lift-curve slope, per radian unless specified 
otherwise 

rolling-moment coefficient, Rolling moment/qSb 
damping-in-roll coefficient, ~ c l-p 

rolling-moment coefficient due to rolling per unit helix 
angle (radians) at tip 

chord, measured parallel to plane of symmetry 
average chord, S/b 
section lift coefficient 

section lift -curve slope, per radian unless specified 
otherwise 

section lift-curve slope in incompressible flow, per 
radian 

distance from leading edge to center of pressure due to 
aileron deflection, fraction of chord 

loading coefficient, cc^/c 

loading coefficient for additional lift distribution 

loading coefficient for basic lift distribution 
loading coefficient for unit-rolling lift distribution 
loading coefficient for residual lift distribution 
loading coefficient for unit angle of attack 
loading coefficient for unit linear antisymmetric twist 
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T]* 


F 

f 

I 


J 

K 


k 


0 


k l 

k 2 

k 3 

k 4 


l 

A 


control deflection in plane perpendicular to hinge line, 
radians 

section lift-curve-slope ratio, c ^ 

dimensionless lateral distance from wing root, 

Lateral distance 

b w 
2 " 2 

plan-form parameter, A/r| cos a 

additional-lift-distribution component-due to sweep 
contribution of function f . to rolling moment due to 
rolling, 4 


abscissa of centroid of area of function f 

correction factor for effectrof taper on lateral center 
of pressure in steady roll 

finite-span correction for wing lift-curve slope, 
equation (7) 

finite-span correction for wing lift-curve slope according 
to slender -wing theory 

finite-span correction for basic lift distribution, 
equation (l2) 

finite-span correction for lift distribution in roll, 
equation (15) 

finite-span correction for residual lift distribution, 
equation (22) 

finite^span correction for rolling -moment coefficient due 
to roll, equation (17) 

section lift 

angle of sweepback at quarter-chord line 


fy* 2 dy* 
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X 

M 

AP 

q 

s 

t 

w 

e 


y* 

y* 



Subscripts: 

II 

III 


effective angle of sweepback in compressible flow, 
tan' 1 ■ ^ A 
- M 2 


taper ratio, Tip chord/Root chord 
free-stream Mach number 

pressure difference between upper and lower surfaces 

dynamic pressure 

wing area 

wing thickness 

width of fuselage 

trigonometric variable corresponding to y*, cos _1 y* 

dimensionless lateral ordinate, Lateral ordinate/Semispan 

dimensionless lateral ordinate of wing center of pressure 

effective lateral center -of-pressure location of resultant 
load causing rolling moment due to rolling 


two-dimensional flow 
three-dimensional flow 


DESCRIPTION OF THE METHOD 
Symmetric Lift Distributions 


The lift distribution for any symmetric angle-of -attack distribu- 
tion may be considered to consist of two parts, a basic lift distribu- 
tion with zero total lift and an additional lift distribution. The basic 
lift distribution for a given twist is defined as the distribution for 
the given wing with the angle of attack reduced equally at every point 
until the total lift is zero. The additional lift distribution is defined 
as the distribution which the wing would carry if it were untwisted and 
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the lift coefficient were equal to 1.0. If the additional and basic 
lift distributions are defined in terms of their loading coefficients, 
then 





and the lift distribution for any variation of angle of attack or twist 
across the span may be written in the form 


y = c i tx a ’ , a + y-b 


(2) 


This section is concerned with the determination of the quantities 
C L 0 ,} 7a? a, and y^. If C^, Cj, and 7 a are known they may, of 

course, be used instead of the values given herein. 

Lift-curve slope.- The wing lift-curve slope may be obtained from 
the section lift-curve slope and a finite -span correction as 


C 


% 


cos A kQ 


( 3 ) 


where c^ and kQ are determined as follows: 

The section lift-curve slope is taken for the section perpendicular 
to the quarter-chard line at a Mach number equaLfo M cos A; it may be 
estimated from the relation 


c l 


C 7 = C7 

L a ‘ai cj 


a 




( 4 ) 
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where the ratio is given in figure 1 from the average of the 

c 7 

theoretical data of reference 6 for airfoils of the NACA 63A, 64 a, and 
65A as well as the 63, 6k, 65, and 66 (with subscripts) series as a 
function of the effective Mach number M cos A for several airfoil 
thickness ratios (perpendicular to the quarter-chord line). For all 
commonly used airfoil sections the lift-curve slope in incompressible 
flow c 7 is known or may be calculated. 

From the value of c 7 a ratio 





(5) 


may be calculated and, hence, a plan-form parameter, defined in refer- 
ence 4 as 


F = 


A 


TJ COS A 


( 6 ) 


According to reference k, the value of kg may then be given in terms 
of this plan-form parameter as 


kg - 



(7) 


The factor kg is plotted in figure 2 as a function of the plan-form 
parameter F. For very large angles of sweepback another factor kg 1 , 
also shown in figure 2, should be used instead of kg. 

Additional lift distribution . - The additional lift distribution 
may be obtained from experimental data or theoretical calculations 
(ref. 7, for instance) or may be estimated from the relation derived in 
appendix A, 



8 


NACA TN 2751 


?a = C l = + C 2 - Vl - y* 2 + C 3 f 
c 


( 8 ) 


where the sweep-correction function f is given in figure 3 and where 
the coefficients Cq, Cp, and C3 are given in figure k. The func- 
tion f depends on an effective angle of sweepback A e defined by 


A 


e 


tan 


-1 


tan A 

yi - m 2 


The elliptic distribution ^Vl - y* 2 is also shown in figure 3 as the 
value of the function f for A e = 0. 

Unlike the value of 7 a given in reference 2, that given by equa- 
tion (8) reduces to the correct value (that given by low-aspect-ratio 
theory or slender -wing theory) when the aspect ratio is very low and to 
the correct (strip-theory) value when the aspect ratio is very high. 

In general, equation (8) may be expected to "apply to most practical plan 
forms with nearly straight quarter-chord lines, except those with very 
large effective angles of sweepback, say greater than about 6o°. For 
such wings the lift distribution may be estimated from the results of 
slender -wing theory (refs. 8 and 9), as explained more fully in 
appendix A. According to slender -wing theory, the lift distribution is" 
likely to be almost-elliptic for a wing with a plan-form parameter less 
than about- 2. 5. A wing with a plan-form parameter greater than about 2.5 
can be divided into three - regions . In the two-dimensional region the 
lift distribution is proportional to the chord, in the tip region the 
lift drops off with infinite slope, and in the root region the lift drops 
to about 2 /n of the value it would have if it were proportional to the 
chord. (See fig. 5.) 

The lateral ordinate of the center of pressure y* of the addi- 
tional lift distribution or of any lift-distribution for a constant 
angle of attack across the span is numerically equal to the moment about 
the origin of the function 7 a , since the area under the function 7 a 
is 1. Therefore, 


y* 


Cl 


1 + 21 
3(1 + 1) 


+ c. 


— + CoJ 

3* 


(9) 
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where J is the abscissa of the centroid of area of the function f 
and is given in figure 6 as a function of the effective angle of sweep- 
hack A e . Equation (9) applies to a linearly tapered wing. F6r a wing 


which is not linearly tapered, the term 
he replaced by 


1 + 2\ 
3(1 + 1) 


of equation (9) may 


_ y*dy* 


Basic lift distribution . - As described in appendix A, the basic 
lift distribution may be obtained from the lift -curve slope and the 
additional lift distribution as follows: 


?b = k l C rJ a - a )7 a 


( 10 ) 


where a is evaluated from the integral 


a 



dy* 


( 11 ) 


and the coefficient k^ is obtained from the relation 


k l 


Ftfl + — + 2 
F 2 


Fl 1 + — + 6 

V F 2 


( 12 ) 


which is plotted in figure 2. The integration indicated in equation (ll) 
may be performed graphically or numerically (by means of Simpson's rule, 
for instance). 
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If there are discontinuities in the angle-of^attack distribution, 
they should be faired before the angle- of -attack distribution is used- 
in equations (10) and ( 11) . Apparently the best results are obtained, 
on the average, when the fairing extends about 0.3 semispan on either 
side of the discontinuity and passes through the midpoint of the dis- 
continuity; the faired curve should have the same area as the unfaired 
one . 


Antisymmetric Lift Distributions 

The lift distribution for_any antisymmetric twist may be considered 
to consist of. two parts: a rolling -type- distribution, which is the dis- 

tribution for the given wing with a linear antisymmetric twist of suffi- 
cient magnitude to have the same rolling moment as the twist distribu- 
tion of interest, and a residual distribution, which is the difference 
between the rolling-type and the true distribution and which, by defini- 
tion, has no rolling moment. If the unit-rolling type of distribution 
is defined by its loading coefficient 




cc l 


where is the rolling-moment coefficient of the distribution for a 

linear antisymmetric twist with unit angle at the tip ^and is, hence, 
the negative of the conventionally defined rolling-moment coefficient 
due to rolling and if the residual distribution is defined by its 

loading coefficient 



then any antisymmetric lift-distribution may be written as 


7 = c Z d <W d + 7 r 


(13) 
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where ae is the tip angle required for a linear antisymmetric distri- 
bution with the same rolling moment as the distribution of interest. 

This section is concerned with the calculation of 7 ^, <x e , 

and 7 r . If the values of 7 ^ and are known they may, of course, 

be used instead of the values presented herein. 

Unit-rolling distribution .- The lift distribution f or a linear 

antisymmetric twist with unit angle at the tip may be obtained from 
reference 10 or 11 or from the relation 


?2 = k 2 C L a ^^a 


(14) 


where 


ko — 


Fll 


+ 4 


+ 2 


(15) 


Fl/1 + ~ + k 


is derived in appendix A and shown in figure 2. The unit -rolling distri- 
bution 7 d is equal to 7 2 divided by the damping- in-roll coeffi- 
cient C^. 

Damping- in -roll coefficient .- The damping- in-roll coefficient may 
be obtained from reference 10 or 11 or, as in reference k, from the 
relation 


C Z d = g- c l a cos A k 4 ( l6 ) 
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where is a factor defined, by 


hi, = 


(17) 


FUl 


l6 . 

+ ~2 + 4 


which is plotted in figure 2. The factor K is a correction for taper 
effects on the lateral center of pressure introduced in reference 12 and 

^L* 


is equal to twice the factor 


b/2 


used in that paper. Another expres- 


sion for — C f may be obtained from the method of the present paper: 

d 


i C7 cos 
4 a 


8 



Co + IC- 


(18) 


where C-^, Cg, and C^ are the factors used in equation (8) and I is 
the moment of inertia of the function f defined as 


I 



fy* 2 dy* 


(19) 


The value of I is given in figure 6. If the wing does not have a 
linear taper the expression 4 f £. jr*2 dy * 

must be substituted -for 


'0 


, . , 2 1 +- 3X 

the term — « 

3 1 + X 


A comparison of equations (16) and (l8) reveals that the factor K 

1 3 1 + X 1 


is equal to the expression 


+ Co H — ICo . 
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^L' 

The two taper corrections K and — are compared in figure j{ a) 


for unswept wings ^since the correction 


b/2 

^L' 


h/2 


does not apply to swept 


wings j and excellent agreement is seen to exist, except for plan-form 
' f L' 


parameters less than 6, for which the values of 


h/2 


given in refer- 


ence 12 are somewhat uncertain. Figure 7(l>) shows the effect of sweep 
on the taper correction factor used in equation (l8). This figure 
serves to explain the statement made in reference 12 that experimental 

^ L P 

evidence indicates that the effect of sweep on the factor — is 


h/2 

likely to be small, a statement which is difficult to reconcile with 
the great deviation of the additional lift distributions of swept 
(particularly sweptforward) wings from elliptical and also of their 
linear antisymmetric lift distributions from that of an. elliptic wing. 
As may be seen from figure 7("b), however, the factor K, which is twice 

and is obtained from spanwise lift distributions, is 


the factor 


b/2 


indeed substantially unaffected by sweep for angles of sweepback between 
0° and about 35 ° , which are the ones for which most of the aforementioned 
experimental evidence was obtained; however, for sweptforward wings and 
for highly sweptback wings, sweep does have the expected effect on K. 

Inasmuch as the correction of equation (l8) permits sweep to be 
taken into account, whereas that of equation (l6) does not, the use of 
equation (l8) appears to be preferable for swept wings. For unswept 
wings both corrections give almost identical results. 

Residual lift distribution .- The residual lift distribution may be 
obtained, as shown in appendix A, from the relation 


7 r = k 3 C j 


■ 


a e) 7 d 


( 20 ) 


where the effective tip angle a e is defined by 
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or 


ot- 



ar d dy* 



7 d y*dy* 

y* 


> 


J 


( 21 ) 


and. where the factor k^ is defined by 



Any discontinuities in the angle -of- at tack distribution must be 
faired before the distribution is used in equation (20) or (21). A 
convenient procedure for fairing this distribution is to plotrthe 
ratio a/y* over the span and to fair it as suggested for discontinuous 
symmetric distributions . The faired distribution may then be used 
directly in the second form of equation ( 2l) . 

ILLUSTRATIVE EXAMPLES 

The method presented in the preceding sections has been applied to 
the calculation of lift distributions for a variety of plan forms in 
order to illustrate certain trends as well_as ..to. compare the results 
with those of other methods. 

Additional lift distributions were calculated for six wings of 
different plan forms in incompressible flow'. The results are shown in 
figure 8. The aspect ratios, angles of-sweepback, taper ratios, and 
plan-form factors of these wings are given in table 1. In the calcula- 
tion of the plan-form factors a section lift-curve slope of 2n was 
assumed. Also given in table 1. are .the constants C d , C 2 , and C^, 

which were taken from figure 4. The lift distributions calculated by 
the method of reference 13 (with 15 points on the span) and, in two 
cases, the lift- distributions calculated by the method of reference l4 
(with 126 lifting points on the span) are also shown In figure 8. There 
is good agreement between the approximate lift distributions given by 
the method of this paper and those calculated by the theoretical methods. 
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The lift -curve slopes and lateral centers of pressure of the addi- 
tional lift distributions of the six wing plan forms considered in fig- 
ure 8 are given in table 1, as calculated from equations (3) and (9) of 
the present paper and as calculated by the method of reference 13 (with 
15 points on the span). The values of kQ required in equation (3) 
for calculating Cq^ were obtained from figure 2 and are also included 

in table 1. Again, there is good agreement of the approximate values 
of Cq^ and 7* obtained by the method of this paper with the theo- 
retical values. 

Also calculated for plan forms 1 and 2 were basic lift distributions 
due to a unit parabolic twist and due to a unit deflection of a half- 
semispan flap (cog6 cos A = l). The results of these calculations are 

given in figures 9 and ID. For the parabolic twist a = y*^, equa- 
tion (ll) yields, in conjunction with the functions / a = taken 

cC L 

from figure 8, a value of a. The function a - a is shown by the 

first curve in figures 9(a) and 10(a). The values of kq were taken 

from figure 2 and are 0.57 and 0.66 for plan forms 1 and 2, respectively. 
These values of a and kq, the approximate values of Cq^ given in 

table 1, and the additional lift distributions given in figure 8 were 
used in equation (10) to calculate the basic lift distributions 

/ cc A 

=(-zr) (.shown by the second curve in figs. 9(a) and 10(a)). 

\ c /b 

The basic distributions given by the second curve of figures 9(a) 
and 10( a) and the additional distributions given in figure 8 may be 
combined with the calculated values of a and Cq^ to give the loading 

for the case where the wing has zero angle of attack at the root and is 
twisted parabolically so as to have a unit angle at the tip. This case 
is illustrated by the last curve of figures 9(a) and 10(a). These par- 
ticular results are compared with the lift distributions which would be 
obtained from the method of reference 13 . The agreement is better in 
this case for plan form 1 than for plan form 2. 

For the effective twist due to flaps, the angle-of-attack distri- 
butions are shown in the first part of figures 9(b) and 10(b). The 
values of a for the faired distribution were calculated from equa- 
tion (ll) and are 0.26 for plan form 1 and 0.27 for plan form 2. The 
curve cc>f a i r ed. - ® was used in equation (10) to calculate' the basic 

lift distributions (shown in the center part of figs. 9(b) and 10(b)) 
and the total lift distributions (shown in the right part of figs. 9(b) 
and 10(b)) as for the parabolic -twist case. The values of Cq^ and kq 
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used to calculate the basic lift due to parabolic twist were also used 
in this case. The lift distribution calculated by the method of refer- 
ence 13 with 15 points on the span (a correction for the. discontinuity 
in the angle-of -attack distribution was included) is compared with the 
lift distribution calculated by the method of this paper. Again, 
better agreement of the approximate with the theoretical distribution 
is obtained for plan form 1 than for plan form 2 . 

In order to illustrate the treatment of antisymmetric angle-of- 

attack distributions, calculations have been made for _ plan form 2 of 
antisymmetric lift distributions due to a linear antisymmetric-twist 
(damping in roll) and due to - a unit deflection of half-semispan out- 
board ailerons ^oq 8 cos A = l). For the linear-twist case the distri- 
bution 72 was calculated from equation (l 4 ), with k^ = 0 .80 obtained 

from figure 2 . The coefficient of damping in roll C7 = 0.33 vas then 

d 

calculated from equation (l8), in which the values k^ = O.63 obtained 

from figure 2 and I = 1.15 obtained from figure 6 were used. The unit- 
rolling distribution was then obtained by dividing t>y this value 

of C i d . The distribution is shown in figure 11 (a). Also shown in fig- 
ure. ll(a) is the unit-rolling distribution calculated by the method of 
reference 13. The coefficient of damping in roll .calculated by that 
method is = Q.38. The unit-rolling distribution and the coefficient 

of damping in roll obtained by the approximate method of the present 
paper and by the theoretical method of refe re nce 13 are in good agreement. 

The function a/y* for the aileron- type angle-of-attack distribu- 
tion is shown at the top of figure 11(b); the fairing extends 0.3 semi- 
span on either side of the discontinuity. The value of a e was obtained 

from the second form of equation (2l), with the approximately calculated .... 
CC Z 

function 7 d s given in figure ll(a), and is <x e = 1 . 12 . The 

FC *d _ - 

value of k^ was obtained from figure 2 and is k^ = 0.68. With these 
values for k3 and ae> the function 73 given in figure 11(a), the 
previously calculated value of C 1 ^, and the faired function a/y* 
given in figure ll(b), the residual lift distribution 7 r = was 

calculated from equation (20 ) and is shown in the middle portion of fig- 
ure ll(b). The lift distribution due to aileron deflection was then 
calculated from equation ( 13 ) and is shown_in the bottom part of fig- 
ure 11 (b). Also given is the lift distribution calculated by the method 
of reference 13 with 15 points on the spaiTand with a correction for the 
discontinuity in the angle of attack.' The _ two lift distributions are in 
good agreement. 
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DISCUSSION 


Inasmuch as experimental determination of pressure distributions 
is a very tedious process, complete information concerning the spanwise 
lift and moment distributions corresponding to all angle-of -attack con- 
ditions is rarely available. Hence, analytical methods for calculating 
these distributions are used almost universally in designing airplanes. 
Most of these methods are theoretical (for instance, for subsonic speeds 
see refs. 9 and 13 to 1 6 ). The reliability of these methods is well 
established; for wings without fuselage, nacelles, tip tanks, or 
external stores, they generally furnish very good approximations to the 
true distributions, provided the angle of attack, airfoil thickness, 
and Mach number are not too large. However, most of them are relatively 
time consuming. In order to overcome this deficiency, calculations for 
many plan forms have been made by some of these theoretical methods 
(refs. 7, 10, and 11, for instance), so that no further calculations 
need be made for the angle-of-attack conditions considered in these cal- 
culations. Another way of avoiding the tedious calculations required 
for the theoretical methods is to use empirical methods such as those 
of references 1 to 5j as well as that of the present paper, which con- 
sists in a combination of these methods. 

Compared with the theoretical methods, an empirical method for 
calculating lift distributions has the disadvantage of being less 
accurate on the average, although the accuracy of the results of the 
empirical method Is often adequate. Another disadvantage of empirical 
methods as compared with theoretical methods is the following: Once 

the accuracy of a theoretical method has been established by comparison 
with known results in a few cases, confidence can be placed in the 
results of this method for widely different cases. When an empirical 
method Is used, the degree of confidence that can be placed In its 
results for cases widely different from those for which its accuracy 
has been established is not nearly so high. As far as the method of 
the present paper is concerned, this limitation implies that for the 
purpose of calculating lift-curve slopes and additional lift distribu- 
tions the method is restricted to plan forms with nearly straight 
quarter-chord lines and effective angles of sweepback which are not 
greater than, say 60°. However, if these aerodynamic characteristics 
have been obtained experimentally or by means of accurate theoretical 
methods, other lift distributions and aerodynamic parameters can be 
calculated by this method for a much wider variety of plan forms by 
using the known lift-curve slope and additional lift distribution; also, 
the accuracy of the final results can be improved in this manner over 
that attainable by starting out with values of the lift-curve slope 
and additional lift distribution obtained by the method of this paper. 
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The matrix scheme of- appendix B and the method for calculating the 
chordwise center of pressure given in appendix C make the method of this 
paper readily applicable to the aeroelastic_ ana lysis of reference 18. 

For wings with straight quarter-chord^lines, moderate angles of 
sweepback, and moderate aspect ratios the lift distributions calculated 
by the method of the present paper have been compared with those cal- 
culated by theoretical methods in connection with the illustrative 
examples. The agreement between the lift distributions and associated 
aerodynamic parameters calculated by the method of the present-paper 
and those calculated by theoretical methods has been noted to be good. 
Inasmuch as the approximate method of- this, paper is based on lifting- 
line reasoning, and inasmuch as lifting-line theory is valid only for 
wings of high aspect ratio, a comparison of the results of the approxi- 
mate method with results obtained by more rigorous theory for wings of 
very low aspect ratio may be of interest. The lift distributions for 
such wings can be obtained from reference 15, provided the wings do not 
have reentrant trailing edges. 

For a wing of very low aspect ratio and with a parabolic twist 

a = y* 2 ; thus = ~^1 - y* 2 , ^ A, and k x = i (see fig. 2); 

also, a, as obtained from equation (ll), is l/4. Hence, from equa- 
tions (10) and (2), 


y-b - f A (y* 2 - 1)4 - y* 2 


and 


7 = |(2y* 2 +-1)^1 - y* 2 
= A^i sin 0 + i sin 3®^ 

(where cos 0 = y*), which is also the result given in table 1 of ref- 
erence 15. That this perfect agreement of the two results is fortui tous 
may be seen from the fact that-for a linear symmetric twist a = |y*|, 
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and 


7 = 





= - A y* + 2- til - y* c 


whereas, according to reference 15, 


7 = 




+ 


sin 


e 


for this case. The two distributions are shown .in figure 12(a) and are 
in good agreement. 

The coefficient of damping in roll of a wing of very low aspect 
ratio can be obtained approximately from equation (l8). From figure 2 
and equation (17) the factor k^ may be seen to be asymptotically 

equal to f /8 as F approaches 0; for F = 0 the factors C]_ and C3 
are 0 and the factor Cg is 1 (see fig. 4), so that as a result of the 
definition of F, 


C 



as given in reference 15. The factor kg = 0.5 for F = 0, so that 


7 2 = Ay*^l - y* 2 

= — sin 26 
2 


as given in reference 15. 
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The -unit-rolling distribution of-a wing of very low aspect ratio 
is then 


7d 



= — y*il - y * 2 


and the factor k^ is l/2. For a parabolic antisymmetric twisty 
a = y~*2 on the pigh-t wing a nd a = -y*^ on the left wing. The effec- 

32 

tive tip angle is a e = * . and the residual lift distribution on the 

15 it 


right wing is 





so that~the lift distribution on the right wing is 


7 



__*2 16 
y* +- 


I5it 




y* 


= A 


8 1 

sin_S + sin 26 + — sin 

15« 8 



whereas, according to reference 15, 


7 


= A 


— cos^g log tanfy- + — \ 

3 it 6 \k 2) 


+ — sin 26 
3 it 


for this case. These two distributions are shown in figure 12(b) and 
are in good agreement. 
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For the cubic antisymmetric twist a = the approximate method 

of the present paper happens to give the same result as reference 15: 


7 = A Q sin 20 + — sin 40 


In general, therefore, the method of the present paper gives results 
for wings of very low aspect ratio which are in excellent agreement with 
the results furnished by low-aspect-ratio theory. 

The method of this paper has certain advantages over theoretical 
methods apart from the fact that the time required for an analysis by 
the empirical method is only a fraction of the time required for an 
a nal ysis by means of one of the theoretical methods and is indeed com- 
parable to the time required to obtain a desired lift distribution by 
interpolating between those furnished in references J , 10, and 11. The 
method of the present paper is rather flexible, so that it may be used 
in some cases where present theoretical methods are inapplicable. For 
instance, in the case of a wing with a fuselage, nacelle, or tip tank, 
neither the method of this paper nor any generally available theoretical 
method can give the lift -curve slope or the additional lift distribution. 
However, if this information has been obtained experimentally, the method 
of this paper can give the lift distribution for any symmetric or anti- 
symmetric twist with an accuracy sufficient for many purposes, (if the 
wing lift -curve slope is different from the lift -curve slope for the wing 
alone, the factors kp, kp, k^, and can be obtained for the value 

of the plan-form parameter which corresponds to the given lift-curve 
slope rather than for the value of the true plan- form parameter. ) 

Another advantage of the empirical method over some theoretical 
methods is the ease with which it lends itself to the calculation of 
aerodynamic influence coefficients, which are required in an aeroelastic 
analysis. (See appendix B.) Also, for a wing in transonic flow the 
theoretical methods used with linearized-theory corrections for com- 
pressibility effects tend to give unsatisfactory results, as pointed 
out in reference 6. The empirical method of the present paper, however, 
is capable of furnishing useful results in certain cases of transonic 
flow because it is based on an interpretation of the three-dimensional 
Glauert-Frandtl rule which permits the use of known section character- 
istics and reconciles that rule with simple sweep theory. (See ref. 4.) 
For unswept wings of low aspect ratio and for sweptback wings which are 
slender in t he sens e of reference 9, that is, wings which have a small 

value of A\jl - M^, equation 3 has been found to yield a good approxi- 
mation to the observed values of the lift-curve slope in some instances 
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of transonic flow up to and "beyond a Mach number of 1, provided the 
wings have thin sections and are at low angles of attack. 


CONCLUDING REMARKS 


A relatively simple method of calculating spanvise lift distribu- 
tions for any angle-of-attack condition has_ been derived by combining 
and extending several previous theoretical methods. The method is such 
that the work can be shortened and the results improved at various stages 
of the computations by introducing experimental or theoretical values of 
certain aerodynamic parameters whenever they are available. In addition, 
the method lends itself readily to formulation of the load-distribution 
problem "In terms of aerodynamic influence coefficients. 

The results obtained by the present method compare favorably with 
those obtained by more time-consuming thedfies. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., May 1, 1952~~ 
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APPENDIX A 

DERIVATION OF THE METHOD 
Additional Lift Distribution 

In the method of reference 1 the additional lift distribution j a 

for unswept wings is given by the average of the wing plan form and an 
ellipse of equal area, so that 


?a = 


cc l 

cC L 



(Al) 


The corresponding relation for swept wings is given in reference 2 as 




7 a f # + f 

2 \ c “o i 


(A2) 


where hIq and hJq are, in the notation of reference 2, the section lift- 
curve slope and the average section lift-curve slope, respectively, and 
where f is a function of y* which is different for each sweep. This 
function f fulfills a purpose similar to that of the function f used 
in the present paper but differs from it numerically. 

That both of these relations are inapplicable to wings of either 
very high or very low aspect ratio may be seen from the fact that for 
very high aspect ratios the additional lift distribution is equal to 
the chord distribution c/c, whereas for very low aspect ratios it is 

defined by the ellipse —ill - y*^ (see refs. 8 and 15), provided the 

3t ' 

wing does not have a reentrant trailing edge. Consequently, to cover 
the entire range of aspect ratios the equation for the additional lift 
distribution should have terms which involve the chord and the ellipse 
and a function which varies with sweep; the extent to which each of these 
three functions enters into the lift distribution should vary with aspect 
ratio. Consequently, the following relation suggests itself: 

7 a = ■C 1 |-+ C 2 ^\|l - X* 2 + C 3 f (A3) 
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where f is a new function which varies with sweep and is given in fig- 
ure 3> and where the C's vary with aspect ratio'^in such a way that Cp 

is 1 for infinite aspect ratio, C 2 is 1 for zero aspect ratio, and the 

sum of all three C's is always 1. A seh_of — C's obtained by analyzing 
the results of the calculations by the method of reference 3 for a large 
variety of plan forms is presented in figure 4. The plan-form parameter 

F = — (proposed in ref. 4) was chosen as the abscissa rather than 

T] cos A ■■■ 

the aspect raticuproper in order to take sweep and compressibility effects 
into account in a manner similar to that employed' in reference 4 for the 
wing lift-curve slope. .... 

In general, equation (A3) may be expected to apply to most practical 
plan forms with nearly straight quarter-chord lines, except those with 
very large effective angles of" sweepback _A e , and to furnish more accurate 

results than do the approximate - formulas of references 1 and 2, inasmuch 
as equation (A3), unlike those formulas, takes aspect ratio into account. 

Wings with very large angles of sweepback can be analyzed by slender- 
wing theory. (See refs. 9 and 1 6 .) The type of wings analyzed in ref- 
erence 9 is characterized by a constant chord which is filleted generously 
at the trailing edge to the extent that~the root chord is jt/2 times the 
chord far away from the root. As may be deduced from reference 9 , this 
wing can be divided into three regions in the manner indicated in fig- 
ure 5. The pressure distribution in the part of the root region forward 
of the apex of the trailing edge is the same as it would be if the wing 
consisted of that part of the root region only, and therefore it can be 
obtained from the delta-wing theory (ref. 8). The spanwise lift distri- 
bution is constant and equal in magnitude to that - in the two-dimensional 
region. In the two-dimensional region the chordwise pressure distribu- 
tion is the same as in two-dimensional -flow; the spanwise lift distri- 
bution is constant - and is given by 


l = qcci a cot A 


In the tip region the pressure is zero behind the leading edge of the 
tip chord, and the spanwise lift distribution drops off to zero in an 
approximately elliptic fashion. For an unfilleted wing the lift near 
the root must be less; if it is assumed to be proportional to the chord, 
the lift at the root is 2/jt times the lift in "the two-dimensional 
region. 

As shown in reference 9; the lift-curve slope depends on the param- 
eter A/cot A, which for large angles of.. sweep is approximately equal 



.M 


naca ffN 




to A/cos A. This; Yal^jc xis $£&&£ lie: appr oMmat ely equal- -tc ‘.-thk- plan'-xo 
form parameter F, "because for a wing with a large angle of sweepback 
the flow in planes perpendicular to the leading edge is almost incom- 
pressible at all subsonic speeds , sd that the section lift-curve slope 
approximately 2 jt. The variation of - C& with A/cot A given in 

figure 7 of reference 9 is therefore reproduced in figure 2 of the 
present paper as a variation of kQ' with F, where kg’ represents 

:: the. ratio: accas .A-’aiior ^lelfdfehi^gT&rf ’fbbn&V rf 91 * snxnoassn 

■ yas tot fioiJifixiJaxb ilxl arid v"icerfi anxi-aniJlil od 

iioidslei erfd yd nsvi^ ax n .^djjdxnd 

Basic Lift Distribution 

In reference l.the basic lift^ distribution 7-5 is given as one- . 
half of the strip-theory ve&Lufe f ~ d&vf - ^ 


' nscra erfd aeomoc xdb / ? £ zqv -Inns dxnxr sjoJ 

^ m \r\ =5 


noxdndx'idaxb orii ,b 


a 


m 


where 

• V l f r\ f 


{t** - x )? „J= ° X< 


a = 


/ ~ a dy* 

'0 c 


(A5) 


- 1 - : -f qoiqqs jfoadde To aelgne beoxrbnx add sue x pD bn 3 jso sno: fw 

In 'rS'f ereh(?d' tr 2 e tSd “statement r is 'niacixF th^t fee wfng ? fffe-ciirve f slope* 

^ ^ /.n b^dmoo^d 


often gives better results in equation (A^) than the section lift-curve 
slope c 1 } so that 


( fOOj-T - ^ = 7 

rb^'C^a-a)! 


(A 6 ) 


Both results were" -Obtslhfed'^if ’h. ‘futely etnp^ic^I^bks T iV/' r 02 

ence 17; however, indicates the possibility of a more rational approach, 
in that it proves by lifting-line theory that an average angle of 
attack a defined by 


*” r b ( f oj-\ 


5 . £l_ 

,-, ci ^ i-, 

x-' w ) -V + "* v A • 


(A7) 


1 
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can be obtained for any symmetric twist distribution from 


a = 


ay a dy* 


(A8) 


The reasoning behind this equation is somewhat as follows: According 

to lifting-line theory the lift distribution for any angle-of -attack dis- 
tribution is given by the relation 



(A9) 


and the distribution for a unit angle of attack across the span by 


n - <=vf(i - -11) 


(A10) 


where and are the induced angles of attack appropriate to 

the lift distributions 7 and 7^, respectively. These equations may 
be combined to yield 


7 - 7i« + (rai x - 7354.) 


so that the lift coefficient for the given twist is 


Cl = f 7 by* 



7pct by* + 




(All) 
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The second integral may he shown to he zero hy setting 



sin nS 


y i = SZ an i 


sin n 9 


(where 9 s cos"ly*) so that, hy lifting-line theory, 


1 r — sin n 6 

“ i = 5V nan ^ 


ai l = d SI 


na 


*1 


sin ng 
Bin 9 


Then hoth the integral of 7<x^ and that of 7j&± are equal to 

na n a ni for all symmetric angle-of- attack distributions, so 


that their difference is zero. Consequently, 


C L = / ?l a d 7* 


(A12) 


and since the additional lift distribution may he obtained hy dividing 7j_ 
hy its lift coefficient (which is Ci^, since a = 1 for the additional 
lift distribution) equation (A8) follows from equat ion (A12) . 

The preceding development suggests an equation for the basic lift 
distribution in the form of equation (A 6), except that the factor c/t? 
must he replaced with cc^cCl in order that the total lift of the 
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bssic iift distribu^^vb^^r^ o^en^e^ . 

tributxon may be written aa ...' 



(A13) 


where kq is a constant which depends on the plan-form parameter F 
and takes the place of the fgctQr, TlJ-2. in equation (a6) . 

i ~ ' ' 

The value of kq has been derived in. reference 3 by means of series 

expansions. In the derivation of reference 3 it- was required thst the 
lifts and rolling monpnt-SuOfA:4ihe-^f-&lIiatribtt-fcifc)n8- given tby <rqi&tion:\(Ai3) 
be equal to the values given by lifting-line theory. In the following 
paragraphs the same result is obtained in a slightly different manner 
which facilitates the extension of t he m ethod to antisymmetric lift- 

distributions. — J -. art 7 — - = ..jo 

B nia n A-* 1 1 

Inasmuch as the values of the lift-curve slope and damping in roll 
as given in references 4- and 12 are based on lifting -line results for 
elliptic wings, a reasonable value for kq should be obtainable through 

analysis of elliptic wlng ) i ! .--?-fer rr ^o i7 r ^i n^s F lifting-line theory pre- 
dicts a lift distribution B riia 1 q M _ £' L 


l iJ r ^ 

7 


-£ 


' n J 

b„ sin n0 

A + 2nr) n 


(A14) 


/Oi-liaiii . AilBtf d£.- fO-sISltS OXXi2Kii~V.3 jX3 UOi 


ci k \L jn9rrpsanoO loias al 9onaT:9"r'ixb xxsri.t 

where q = — — and b n is defined by 
2m 


(SIA) 



sin n0 


(A15) 


r ; 


% 


gnxhlvxb yd benxsido 9d ysm noxiudliiaxh_ fill Xsnoxiibhs 9hl aorua a- 
enoxilbba 9dd dot i . = .,£> aooia Aj> ip ai rioiriw) JnaioiTiaoo i'll I sj; 


► (SIA) noldsnp9 moii awciiol t&A) noiiBupe (.noldird trial*. .ft i 

■fill otsad edi poi rferi ??«Tp 9 / n<5 ct 3 in ; aiSqjJ&y si yarbaoarrq 9 . IT - 
TT \ 0 lojosl 9dJ derii dqsoxa (,dA) noiisnpe To nrro'i srid ni noj-Jj/o i 
*><$ to d’ixl Isiod arid Jsdi tstno- “ni i^-tO o ^ j oo rfiiw baoslqoi ■ ■ 
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For a symmetric angle-of -attack distribution, only odd values of n 
are needed in equation (Al4) . The lift coefficient is then given by 


C L » | =1 - an, b x (A 16 ) 

Since the lift is to be zero for the basic lift distribution, the first 
term, of the series must be omitted, so that for the symmetric case the 
first term of the series is that for which n = 3- For the case of 
constant angle of attack, however, b^ = a and all other b's are 
zero, so that 


8nA . 

7 - 1 = * a sin 6 

1 A + 2t] 


(A17) 


'Itt 


A + 2rj a 


(A18) 


and 


4 

■st 


7 a = — sin 6 


(A19) 


Consequently, substituting equations (A15), (Al8), and (A19) into 
equation (A13) gives 


7b = k l 


A + 2q ° 1 cl it 


b n sin n0 


(A20) 


since a = b-j_. This distribution should be the same as that given by 

equation (Al4) (starting with the term n = 3). Clearly, no value 
of kj_ can satisfy this condition exactly. If the first terms in the 

two series are set equal, however, the others will not differ much; 
furthermore, since the first term is usually much larger than the others, 
the lift distributions will generally be almost identical if the first 
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terms of the series are equal. In order to satisfy that condition, 
k]_ must be taken as 


A + 2t| 
A + 6tj 


(A21) 


which is identical to the result obtained in. reference 3, except that 
there q was assumed to be 1. - — 

A lifting-surface and sweep correction for k-j_ may be deduced in 

the manner employed for the damping-in-roll .coefficient-on references 4 
and 12. * If the product k]_C^ that occurs in equation (A13) is obtained 

from equations (Al8) and (A2l) and written in the form 


k l C Ia- = a 


A/3 


c l 


f + 2n 


a 


a correction identical to that for Cj^ but based on one-third the 

aspect-ratio suggests itself. Consequently, upon introducing the sweep 
and lifting- surf ace corrections in this equation and dividing by Cj^, 
the following relation is obtained: 


k l " 


\ 1 + u(3_S| LA) 2 + a, cos A 


(A22) 


A^l + 36 ^— CQS — ^ + 6t) cos A 


This factor is plotted in figure 2 against the plan-form parameter 
A 


F s 


T) COS A 


Lift-Distribution in Roll 

If an equation similar to equation (A13) is applied to a linear 
antisymmetric angle-of-attack distribution" with unit value at the 
tip (a = y*), the relation 
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7 2 = k 2 C i^y*7 a (A23) 

is obtained. A factor k 2 is 1136(1 rather than k]_ because of the dif- 
ferent nature of the induction effects for the symmetric and antisym- 
metric cases. , The desired value of k2 may 1)6 obtained from the lift 
distribution of elliptic wings , as follows: 

For the angle-of-attack distribution in this case , equation (A15) 
yields b 2 = i, and all other b's are zero. Hence, equation (Al4) 
yields 




^ i sin 26 


A + ^1 2 


whereas equation (A23) yields, after substitution of equations (Al8) 
and (A19), 


7 2 = * 2 c l 


a A + 2 t] jt 


— sin 9 cos 6 


The desired values of k 2 
sions for 7 2 : 


\ 

^ay be obtained by equating these two expres- 


k 2 


A + 2tj 
A + 4t) 


(A24) 


This value, also, is identical to that obtained in reference 3 for tj = 1. 
It may be corrected for lifting-surface and sweep effects in the same 
manner as k-^, so that 


k 2 = 


AU1 + 4j 


r 

fr\ cos A y 

\ A / 


+ 2t\ cos A 


A| 1 + j + 4rj cos A 


(A25) 


This factor is also plotted in figure 2 against the plan-form parameter F 
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Damping- in-JRoll Coefficient 

By calculating the rolling moment of the lift distribution given by 
equation (A23), a value for the damping-in-roll coefficient ~ c Zp' 

may be obtained: 


° l d " | 


‘-a 


cos A 



7 a y* 2 dy* 


(A2 6 ) 


The product of the factors kg and kQ is referred to aB the factor k^_ 
in the following derivation. 


If equation (A3) is substituted into equation (A26 ) , the damping 
coefficient may be written as 


C l d ’ 5 ““ A **(! Cl + C2 + IC 3) 


(A27) 


where \ is the taper ratio , where I is_ defined by 


1=4/ fy* 2 dy* 
'■'0 


(A28) 


and is plotted in figure 6 (b), and where the expression 


4 



— y* 2 dy^ 
c 


must be-substituted for the term 


2 1 + 

3 1+ J. 


for wings which do not~have 


linear taper. The value of Cj given by equation (A27) contains a 
correction for taper (the Cj_ term) and a correction for the effect 
of sweep on the spanwise lift distribution (the C 3 term). 
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Antisymmetric Lift Distributions 

Inasmuch as equation (A23) has been derived specifically for linear 
antisymmetric angle-of-attack distributions., it may not yield accurate 
results when applied to other antisymmetric angle-of-attack distribu- 
tions, as is done in reference 3- The derivation of the equations for 
symmetric lift distributions, in this appendix suggests a parallel deriva- 
tion for antisymmetric distributions which has as its starting point 
not the relation given in reference 17 , but rather its analog for anti- 
symmetric lift distributions. The resulting equations for antisymmetric 
lift distributions are not as simple as equation (A23) - they imply a 
separation of the lift distribution into a damping-in-roll and a 
residual lift distribution analogous to the additional and basic lift 
distributions in the symmetric case - but they are applicable to a 
wider variety of antisymmetric angle-of-attack distributions than is 
the method of reference 3 • 

According to lifting- line theory applied to the linear antisymmetric 
case (a = y*), 

r 2 * c laT( y * • “is) (A29) 


where a± 2 is the induced angle of attack corresponding to the distri- 
bution 72 « Combining this equation with equation (A9) for any arbitrary 
lift distribution yields the relation 


77 * = 72 a + (7ai 2 - 72 a i) 


(A30) 


The rolling-moment coefficient of any antisymmetric lift distribution 7 
may be written as 


C 


l 



7y*&y* 


so that, from equation (A30), 


(A31) 



'0 


7pa dy* + i J ( 7 a i2 - 7 2 a i) <ly* 
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By reasoning similar to that employed in connection with equation (All) 
the second integral may be shown to be equal to zero for any antisym- 
metric 7 distribution, so that 





7 2 a 


c 


z 



a 




(A 32 ) 


A damping- in-roll distribution 7^ for unit rolling- moment coef- 
ficient may be defined in a manner analogous to the additional lift dis- 
tribution as 


?d = 


l2_ 

C *d 


since Cz d is the rolling-moment coefficient of the 72 distribution. 

An effective angle of attack 0^, which is the angle of attack at the 
tip of a linear antisymmetric distribution with the same rolling moment 
as that of the given antisymmetric distribution 7 , may then be defined 
as 


or 




7 d a dy* 



(A33) 
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The analog of the basic lift distribution (which has zero lift) in 
the case of an antisymmetric distribution is the residual lift distri- 
bution 7 r , which has zero rolling moment. Equations (A13) and (A32) 
then suggest the form for this distribution: 


?r 



“p \7 f 


(A3b) 


The rolling moment of y v is then zero, as required. 


The factor may again be evaluated through an analysis of 

elliptic wings by lifting-line theory. For the antisymmetric case, only 
even values of n are contained in the series of equation (Al^) . Fur- 
thermore, since 




■A 
A + 



‘ (A35) 


then for y r the term with n = 2 must be zero, and the series must 
begin with the term n = 4. But for an elliptic wing 


1 A 

C ^d 8 Cl co A + 4 t] 


(A3 6) 


?d 


— sin 26 


(A37) 


and 


CL. 


= 2b r 
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so that- equation (A3*0 becomes 


7 r = ko ~ ci - . 4 - Y — - 2b^\ — 
r 3 8 ia A + J^yy* it 


sin 20 


or, as a result^ of equation (A15), 


= k 3 8 C >a 


8 l a A + 4q 


— ^ b n sin n0 

11 4 


(A38) 


The values of y v given by equations (A38) and (Al^) cannot~be exactly 

equal for all antisymmetric angle-of-attack conditions, but if the first 
term of the series (which is usually much larger than all others) is made 
the same for both, the others will be approximately equal. Hence k^ 

must have the value 


k, . t±±L. 

J A + 8t) 


(A39) 


which, by the same reasoning as that for k]_, may be corrected for 
lifting-surface, sweep, and compressibility effects by setting 


cos A 


k 3 = 


4 + - M 

l^-_ 64 ^— .- os - — ^ + 8r) cos A 


(A4o) 


This factor is plotted in figure 2 against the plan-form parameter F. 
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APPENDIX B 

CALCULATION OF AERODYNAMIC INFLUENCE COEFFICIENTS 


Aerodynamic influence coefficients are useful in methods of static 
aeroelastic analysis, particularly those employing matrices, such as 
that of reference 18. These influence coefficients may be considered 
to represent the lift (in dimensionless form) at one point of the span 
due to a unit angle of attack at another point. Such influence coef- 
ficients can be obtained by means of some of the theoretical methods 
by calculating the lift distributions corresponding to angle-of-attack 
distributions which are zero everywhere on the span except for one point, 
at which they are infinite in such a way that the area under the angle- 
of-attack distribution is zero. However, such lift distributions are 
difficult to calculate, and they may be obtained at points along the 
span which are inconvenient for an aeroelastic analysis. 

The influence coefficients obtained by the procedure of reference 5 
are based on the method of calculating spanwise lift distribution pre- 
sented in reference 2. Individually these empirical coefficients do 
not have any significance, as do the coefficients described in the pre- 
ceding paragraph, but when multiplied by the values of the angle of 
attack at several points on the span and summed, they do yield the lift 
at a point on the span in the same manner as the other influence coef- 
ficients. Apart from the fact that they can be calculated more simply 
and for any points on the span for which a set of integrating factors 
can be calculated, the empirical influence coefficients have the same 
advantages and disadvantages in comparison with the theoretical coef- 
ficients as do empirical methods in comparison with theoretical methods. 

In this appendix a procedure is outlined for calculating influence 
coefficients In a manner which Is similar to that employed in refer- 
ence 5> but which is based on the method of calculating lift distribu- 
tions given in the present paper rather than that of reference 2. In 
essence the procedure presented in this appendix consists in formulating 
this method of calculating lift distributions in matrix notation by using 
integrating matrices of the type presented in reference 18. Compared 
with influence coefficients calculated by the procedure of reference 
the ones calculated by the procedure presented in this appendix have the 
advantages of greater accuracy and wider applicability, because the 
method of calculating spanwise lift distributions presented herein, on 
which the coefficients are based, has these advantages over earlier 
methods. The influence coefficients given by the procedure of refer- 
ence 5 are applicable only to symmetric loadings and wings of moderate 
aspect ratio, whereas those given by the procedure of the present paper 
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ere not subject to either of these restrictions. They should, there- 
fore, be useful for most purposes for which aerodynamic influence coef- 
ficients are required. 

The matrix notation used herein is the. same as that used in 
reference 18 . 


Symmetric Case 

Equations (2) and (10) may be combined to yield 


1 - k lWa + k l a7 a 

or, in matrix notation. 



where is a column matrix which consists of the values of y at 

several stations on the wing, is a column matrix which consists 

of the values of a at those stations, £cify is a column matrix of 
elements all equal to a, and l/al is a diagonal matrix the nonzero 
elements of~which are the values of y a at the same stations. 

The integration indicated by equation (ll) may be written in 
matrix form as 



a Ha root + b fljWfl 


(B2) 


where H is defined by 


H = 


,w/b 


0 


7 a <ly* 


(B3) 
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and may be evaluated with sufficient accuracy for the present purpose 
from the relation 



and where [I^J is an integrating matrix. As described in reference 18, 

a suitable integrating matrix may be obtained by approximating the inte- 
grand (the function a 7 a ) by parabolic segments. In view of the fact 
that 7 a goes to zero with infinite slope a special type of parabola 
has to be used, for instance, one which consists of a linear combina- 
tion of the l/2 power and 3/2 power of the distance from the tip. With 
this approximation to the curve of ay a , the integrating matrix for 
stations at the wing root and every sixth of the distance from the wing 
root to the wing tip, for instance, may be written as 


[ii| = [0.05556, 0.20833, 0.1527Q, 0.16667, 0.14913, 0.22500, o| (b 4) 

Root Tip 


Equation (B2) can be written as 


3 * LhJWp} 


where the row matrix [^IjJ [j^ is obtained by calculating the row 
matrix — ^ - [j-ljl/aj a 11 ' 3 - adding H to its first element. A square 


matrix Ql£] [7^] can then be constructed which consists of rows all 

equal to |I]J [7^] . With this square matrix, equations (Bl) and (B2) 
can be combined and written as 


pa} - Cl^al []l * *l)[ll] W + k lL 1 l]{“} 
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ko 


or 


{7a} = C Ia [Q s ]{a} (B5) 


where [^1~] is the unit matrix., and where .the matrix QQ s } defined hy 


M ■ N R 1 - k i) + k iL 1 _ f) 


(b6) 


is, in effect) an aerodynamic -inf luence -coefficient matrix for symmetric 
lift distributions. - 

In these matrices the values of 7 , y 8 , and a are all taken at 
the same set of stations, and the integrating matrix must be set up for 
the same stations. Stations at the wing root and every sixth of the 
distance from the wing root to the wing tip (but - excluding the wing tip 
proper, where y 8 is zero unless tip tanks or end plates are present) 
have been found convenient for many purposes. For these stations 


t]* = 0, 0.1667, 0.3333, 0.5000, 0.6667, 0.8333 


where q* is the ratio of the lateral distance of the given station 

t) - w 

from the wing root to the length — - — . The integrating matrix given 

by equation (B4) may be used for these stations. The last element, 0, 
of the matrix of equation (b 4) can be disregarded. 


Antisymmetric Case 

The antisymmetric lift distribution may be obtained from a combina- 
tion of equations (13), (20), and (2l) by use of matrix methods in the 
same manner as in the symmetric case. The resultris 


{7} = c *dN 




(B7) 



SM 
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hi 

or 


&} C idpa]“L a 3’ 


(B 8 ) 


where 


M - L’-il 




(B9) 


is 


an 


distributions. The 


“atrix for 


aer ° aroa " C -inf 1Ue nce-coemcient 

equal to - z b 1 r- 1 1 .. -1 ' ~ 1 square matrix with rows all 


entisymmetri 


c lift 


Mtm **a 

b L*ijL«n and with a constant w* 
element. This constant - - - t0 the flrst 


is defined by 


H- . i £ r 

2 V Jn 


w/b 


y*?& dy* 


and is given 


approximately by the relation 


S' 



The matrix 
relation 


!y* 


can be calculated at -hh* „• 

given stations from the 


y* = 



w 

b 



k2 
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APPENDIX C ■■■■ 
MOMENT DISTRIBUTIONS 


The s panwise distribution of the pitching moment is often of less 
interest than the lift distribution; for aircraft loads calculations the 
lift distribution usually is of primary interest^ and for stability 
calculations the use of two-dimensional centers of pressure often yields 
sufficiently accurate over-all aerodynamic parameters. For some pur- 
poses, however, such as static aeroelastic calculations, a knowledge of 
the spanwise distribution of the pitching moment is important^ - ' In this 
appendix some available knowledge concerning this distribution is 
summarized for use in conjunction with the method of this paper. 


The variation of the pitching moment along the span may be obtained 
from the lift, distribution and the local centers of pressure. The local 
aerodynamic center, that is, the center of pressure due to angle of 
attack or twist of an uncambered wing, and the center of pressure due 
to control deflection are of most general interest.' 

Compared with the amount^of information available on spanwise lift 
distributions, relatively little is known about the local centers of 
pressure of the lift distribution. The commonly used Weissinger method 
(ref. 13 ) for calculating lift distributions is characteristically inca- 
pable of furnishing local centers of pressure. Methods suitable for 
obtaining such information, such as Falkner's (ref. 14), are very 
cumbersome, and the local aerodynamic centers calculated by these 
methods are not altogether reliable. For_ instance , by using different" 
vortex representations in Falkner's method different“local aerodynamic 
centers are obtained, although the lift distributions are nearly the 
same. The centers of pressure due to aileron deflection cannot be 
calculated accurately by Falkner's method. An additional uncertainty 
in connection with the local centers of pressure is the fact that these 
centers are more sensitive than the lift distribution to deviations from 
the potential-flow conditions (as a result_of boundary-layer separation, 
for instance) assumed in the analytical methods. 

As a result of these considerations the information concerning 
local centers of pressure presented in subsequent sections should be 
regarded as qualitative rather than quantitative in nature; it is 
intended as a rough guide until more refined methods of calculating 
such information are available. The values of the local aerodynamic 
centers and centers of pressure due to control deflection are for 
incompressible flow. In order that subsonic compressibility effects 
may be taken into account approximately, the values should be estimated 

for_an equivalent wing with an aspect ratio equal to /l - times 
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the actual aspect ratio and with an effective sweep angle A g , the 

tangent of which is l/l/l - M 2 times the tangent of the actual sweep 
angle A . 


Local Aerodynamic Center 

The local aerodynamic centers of an unswept wing of very high 
aspect ratio are at the quarter-chord line. As the aspect ratio 
decreases the local aerodynamic centers move forward, in particular 
near the tip of the wing. For constant -chord wings of vanishingly low 
aspect ratio, virtual-mass considerations indicate that the local cen- 
ters of pressure are near the leading edge, although the applicability 
of low-aspect-ratio theory to such a wing may be questioned. 

The local aerodynamic centers of several unswept constant -chord 
wings calculated by Falkner’s method and those of two elliptic wings 
calculated from the results of reference 19 are shown in figure 13 (a) . 

At high and medium aspect ratios the wing taper has little effect on 
the local aerodynamic center, as may be seen, for instance, by comparing 
the results for the constant -chord and elliptic wings. At very low 
aspect ratios, on the other hand, virtual-mass considerations indicate 
that the effect of taper is likely to be more pronounced. If the leading 
edge is swept back, the local aerodynamic centers of a wing with vanish- 
ingly small aspect ratio may be estimated from the data shown for the low- 
aspect-ratio delta wing in figure 13 (c) by considering the part of the 
wing enclosed by the leading edge and the line which connects the two 
ends of the leading edge to be a delta wing (see fig. 5) • The 
aerodynamic -center locations shown in figure l4(a) were obtained in this 
manner . 

The local aerodynamic centers of swept wings of very high aspect 
ratio are also on the quarter-chord line, except at the root and the 
tip. Near the wing root the local aerodynamic center, which is at the 
quarter-chord point for an unswept wing, moves back as the wing is swept 
back and approaches the midchord position for very high sweep angles; 
as the wing is swept forward the local aerodynamic center moves forward. 

At the wing tip, on the other hand, the local aerodynamic center moves 
forward toward the leading edge as the wing is swept back and rearward 
as the wing is swept forward. 

The local aerodynamic centers for swept wings of medium aspect 
ratio depend both on the aspect ratio (in much the same way as for 
unswept wings) and on the angle of sweep (in a manner similar to that 
described in the preceding paragraph) . The results of calculations by 
means of Falkner’s method for constant -chord wings of various sweep 
angles and aspect ratios, as well as for one tapered wing and two delta 
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wings, are shown in figures 13(b) and. 13(c). Comparison of the two plan 
forms of aspect ratio 2 with 60° sweepback corroborates the previously 
made statement that taper has an effect on the local aerodynamic center 
at low aspect ratios, although the effect-is fairly small. 

For wings of very low aspect ratio with sweptback leading edges 
and without reentrant trailing edges, the method of reference 3 may be 
used. In the case of delta wings it yields the aerodynamic -center 
variation along the span shown in figure 13 (c). 

For wings with angles of sweepback larger than about 60° the 
results of slender-wing theory presented in reference 9 may be used as 
a guide. The local aerodynamic -center positions obtained in this manner 
are shown in figure 5* These variations are only approximate, particu- 
larly in the case of the wings with small plan-form parameter, for which 
only the values at the root and tip are known. 

The local -aerodynamic -center lines shown in figure 13 have been 
calculated for constant angles of attack along the span. From the data 
of reference 19, however, a local-aerodynamic -center line can be calcu- 
lated for a linear antisymmetric angle -of -attack distribution of an 
elliptic wing of aspect-ratio 6. This line agrees perfectly with the 
one shown in figure 13 . Consequently, it may be assumed that the nature 
of the angle -of -attack distribution along the span has little effect on 
the local -aerodynamic-center location for unswept wings, and probably 
for swept wings as well. 


Local Center of Pressure Due ta^Control Deflection 

For unswept wings of high aspect ratio and, except at the root and 
the tips, for swept wings of high aspect _ ratio, the two-dimensional 
value of the center of pressure due to control deflection may be used. 
For medium aspect ratios no results comparable to those for the local 
aerodynamic center are available. In order that some indication of the 
probable location of the center of pressure relative to the two- 
dimensional value may be obtained, the difference between the lift 
actually carried at~any section of the wing and that which would be 
carried. if the section were in two-dimensional flow may be considered 
to act at the aerodynamic center. This assumption leads to the equation 

Acp 6 = cp Sm - cp 5lI 



7ll “ ’'III 


III 
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where 7 m is the value of 7 calculated by the method of this paper , 
a is the local aerodynamic center, and 


= 2. cos A(S cos A) 

c 


c o 

= — c 7 6 cos^A 


where 5 is the control deflection measured in a plane perpendicular 
to t he hinge line and ag is the two-dimensional value of the control 

effectiveness . 

For wings of very low aspect ratio the method of reference 8 may 
sometimes he used to calculate centers of pressure due to control deflec- 
tion, provided both the wing and the control surface have sweptback 
leading edges, the trailing edge is not reentrant, and the hinge line is 
perpendicular to the free stream. The centers of pressure due to a 
symmetric control deflection as given by the apparent-mass concept on 
which reference 8 is based are shown in figure l4(b) for two such wings. 
When the entire tip of a low-aspect-ratio wing is rotated about a hinge 
line perpendicular to the free stream (as shown in the third sketch of 
fig. 14(b)) the resulting spanwise lift distribution can be obtained 
from reference 15 , but no chordwise distributions are given by this method. 

In the case of the first two plan forms shown in figure 14(b), all 
the lift carried ahead of the hinge line when the whole wing is at an 
angle of attack is carried at the hinge line when only the control sur- 
face is deflected. The control surface itself carries the same lift as 
it would if the whole wing were deflected. 

For wings of low aspect ratio the apparent-mass concept also indi- 
cates that the two-dimensional value of 05 cannot be used; instead a 
value close to 1 should be used for wings of very low aspect ratio such 
as those shown in figure l4. 
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TABLE 1.- GEOMETRIC CHARACTERISTICS, CONSTANTS , AND LIFT -CURVE 
SLOPES OP THE SIX PLAN FORMS CONSIEERED IN FIGURE 8 


Plan A A, 
form deg 





c l^ 

y* 

C 2 

C 3 

k o 

Equation 3 

Method of 
refer- 
ence 13 

Equation 9 

Method of 
refer- 
ence 13 

0.44 1 

0.254 

0.721 

4.529 

4.321 

0.430 

0.425 

.266 

.327 

, .792 

3.517 

3.444 

. 44 a 

. .455 

.266 

.327 

.792 

3.517 

3.237 

.471 

.480 

.652 

.163 

.577 

3.140 

2.795 

.449 

.452 

.584 

.194 

.634 

2.817 

2.843 

.4i4 

.407 

.441 

.254 

.949 

2.265 

2.290 

.445 

.455 


NACA^ 

































0 2 4 6 8 10 12 14 16 


Plan-form parameter, F 

Figure 2.- Factors Fq, and 1%. 
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Effective angle of cweepback^A. ej deg 

(a) Spanwise center of pressure J. 



Effective angle of sweepback,Jl e; deg 


(b) Moment of inertia I. 

Figure 6.- Centroid of-area and moment of inertia of lift-distribution 

function f. 
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Approximate method — Method of reference !3 




(b) Unit deflection of half- semi span flap. 

Figure 9 .- Symmetric lift distributions for plan form 2. VJ 
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Approximate method Method of reference 13 



(a) Unit parabolic twist. 



(b) Unit deflection of half- semi span flap. 

Figure 10.- Symmetric lift distributions for plan form 1. 
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Dimensionless lateral ord/nofe, y* 


(a) Linear twist ( damping- in -roll case). 

Figure 11.- Antisymmetric lift distribution for plan form 2. 
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QTd ) Unit deflection of half- semispan flap. 
Figure 11.- Concluded. 
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Approximate method 

-Method of reference !5 


.? .3 A .5 .6 .7 .8 

Dimensionless, lateral ordinate, y* 

(a) Symmetric linear twist. 



0 J .2 .3 .4 . 5 .6 .7 .8 .9 

Dimensionless lateral or d mate, y* 

(b) Antisymmetric quadratic twist. 

Figure 12 .- Lift distributions for wings of very low aspect ratio 
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(a) Aerodynamic- center locations. 
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(b) Center-of -pressure locations due to control deflection. 

Figure 14.- Aerodynamic-center locations and local centers of pressure ctn 

due to control deflection on lov -aspect-ratio vings without reentrant OJ 

trailing edges according to low-aspect-ratio theory. 


HACA TN 2751 





